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I Abstract. I discuss the high resolution XMM-Newton Reflection Grating Spectrometer (RGS) 

spectrum of SN 1006. SN 1006 is one of the best examples of a supernova remnant that is far 
out of ionization equilibrium. Moreover, optical, UV and X-ray data indicate that it is also out 
O i of temperature equilibrium. I discuss the X-ray evidence for this. 

^ • In addition I discuss the lower resolution RGS spectrum of the eastern rim of SN 1006. Despite 

' the lower resolution, the spectrum contains significant evidence for an asymmetric expansion veloc- 

ity. Two likely solutions fit the O VII triplet. One with no significant thermal broadening and a shell 
velocity of ^ 6500 km s^ ' , and one with significant broadening and a shell velocity of 9500 km s^ ' . 
The first solution seems the most plausible, as it is consistent with radio expansion measurements, 
which suggest a decelerated shell. 

> 

^ : INTRODUCTION 
(N 

^ I The study of supernova remnants lies at the crossroads of astrophysics, combining the 

O ■ interests in the end products of stellar evolution with the investigation of the interstellar 

r-| ! medium and cosmic ray acceleration. This is all the more interesting now that Type la 

Qh| supernovae [e.g. 1] have become central to the study of the geometry of the universe and 

Q ■ gamma ray bursts are known to originate from a subset of core collapse supernovae [2]. 

However, of more interest for these proceedings is that supernova remnants are objects 
in which interesting plasma processes take place. Supernova remnant shocks move 
through a medium with a typical density of 1 cm^. At those low densities Coulomb 
(particle-particle) interactions are rare, leading to typical collision time scales of several 
r> I thousand years. Shock heating does, however, occur in supernova remnants, otherwise 

• we would not observe X-ray emission from hot plasmas. This implies that the heating 

process takes place through collective effects, involving plasma waves. Apart from being 
coUisionless, supernova remnant shocks are also the most probable source of cosmic 
rays, at least for energies up to ~ 10^^ eV, but possibly even up to ~ 10^^ eV [6, 7]. 

The fact that supernova remnant shocks are coUisionless and are the locations of 
particle acceleration has two important consequences. First of all, we can no longer 
assume that different particles species (protons and other ions and electrons) are in 
temperature equilibration. In an extreme case, and ignoring particle acceleration, the 
temperature, kTi, of each plasma component, i, is: 



where y is the adiabatic index, m,, is the particle mass, and v., is the shock velocity. 
For full equilibration this is fcr = 3/16 < m > v^. In addition, we have to consider that 




□ectron Temperature (eV) 

SN 1006 Tycho Cas A 

T 1 — I I I 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I I I M| 1 1 — I I I 1 1 1 1 1 1 



He-like Fully ionized 




FIGURE 1. An illustration of the effects non-equilibration ionization. Although the top and bottom 
panels look very similar, the top panel shows the oxygen ionization fraction as a function of electron 
temperature for the case of collisional ionization equilibrium (CIE), whereas the bottom panel shows how 
in the case of recently heated plasmas (NEI) the ionization fraction at a fixed temperature of kT^= 1 .5 keV 
changes as a function of ionization age, which is the product of electron density and time since the 
plasma was heated {n^t). Note the sturdiness of the He-like ion, which is present over a wide temperature 
or list range. The calculations are based on the simple to use, but somewhat outdated ionization and 
recombination rates in [3]. 

efficient cosmic ray acceleration, and possible magnetic field generation, will take away 
energy, lowering the plasma temperature(s) of the thermal plasma components [8]. 

As a result measuring the plasma temperature through X-ray spectroscopy, which 
usually will give us the electron temperature, will no longer inform us directly about the 
shock velocity. This has been known for quite some time [e.g. 9, 10], but until recently 
it was ignored consciously [e.g. 11] or subconsciously, as it was difficult to assess the 
amount of temperature equilibration from the observational data. 

However, another form of non-equilibration, namely non-equilibration of ionization 
(NEI), has received more attention over the last three decades, because its consequences 
could be more easily deduced from X-ray spectra of sufficient quality [e.g 12, 13, 14, 
11]. 

The concept of NEI is relatively simple [9, 15, 16]. NEI is important in supemovarem- 
nants as in the relatively recently heated plasmas the ions have not had sufficient time to 
ionize up to the ionization level at which the number of ionizations is compensated by 
the number of recombinations (i.e. collisional ionization equilibration or CIE). The ef- 
fect of NEI is illustrated in Fig. 1 . Observationally NEI gives rise to a mismatch between 
the electron temperature derived from line emission ratios, and the electron temperature 




FIGURE 2. An illustration of the effect of non-equilibration of temperatures, taking into account 
ionization effects, but not adiabatic expansion. The temperature of the electrons (dotted), protons (solid), 
helium (dashed-dotted) and oxygen ions (dashed) is proportional to their mass ratio behind the shock. 
Due to Coulomb interactions the temperatures slowly equilibrate. It is interesting to note that the oxygen- 
proton equilibration will eventually be faster than the helium-proton equilibration, as the oxygen ions 
acquires a higher charge and the interactions scale with 7} [4, 5]. 



derived from the continuum shape, which reflects the actual electron temperature. In 
addition, spectra of NEI plasmas will display lines that are unique for NEI, and are the 
result of inner shell excitations and ionizations. The Fe-K line emission is most widely 
used as a diagnostic tool for NEI conditions, as for a plasma with 2 keV the cen- 
troid of the n = 2 — 1 line emission shifts gradually from ~ 6.4 keV to ~ 6.7 keV for 
FeXXV [17, 18]. 

Here I discuss high resolution spectra of a supernova remnant that is one of the most 
extreme cases of non-equilibration of ionization and temperatures, SN 1006. 



THE ZMM-iV£WrOiV-RGS SPECTRA OF SN 1006 

The young historical remnant SN 1006 is thought to be a Type la supernova remnant on 
account of its high Galactic latitude (-1-14.6) and the long duration of the historical light 
curve [20]. The high latitude also explains the relatively low density of the local medium 
of 0. 1 cm-3 [21, 22, 23, 24]. This makes SN 1006 the ideal candidate for studying the 
two non-equilibration effects, as it has an ionization parameter of ne?~ 2 x 10^ cm~^s. 
The distance to SN 1006 is 2.1 kpc [25]. 

XMM-Newton observed SN 1006 several times, both as a guaranteed time and guest 
observer target. I will discuss here the guest observer observations of August 2002, 
which consisted of two pointings. One pointing was on the bright knot at the northwest- 
em edge of the remnant, the other was on a narrow feature on the eastern rim (Fig. 3). 
SN 1006 is a large remnant (30'), but the pointings were optimized in order to obtain 
high resolution spectra of bright, small structures. Moreover, the pointings were on the 





FIGURE 3. Chandra mosaic of SN 1006. The circles indicate the XMM-Newton pointings discussed 
here. The emissivity profile is taking in the dispersion direction of the RGS [from 19]. 



rest of the remnant where it is unlikely that there are bulk plasma motions in our line of 
sight. This is necessary in order to allow estimates of the thermal line broadening. 

A sufficiently high resolution spectrum was obtained for the northwestern pointing, 
but for the eastern rim, only the non-thermal X-ray emission is very narrow, whereas 
the thermal emission is much more extended, degrading the spectral resolution. This 
is somewhat unfortunate, as the electron and ion temperatures of the plasma may be 
influenced by efficient cosmic ray acceleration, and non-thermal emission from the rims 
of SN 1006 indicate a very efficient electron cosmic ray acceleration at the west and 
eastern rims, but not at the northwestern edge [26, 27]. 



The X-ray spectrum of the northwestern knot 

The XMM-Newton Reflection Grating Spectrometer [RGS, 28] spectrum of the north- 
western knot shows many spectral details that cannot be seen in CCD spectra. The most 
dominant emission is coming from O VII, in particular from the O VII triplet, but also 
O VII line emission from higher quantum levels up to the series limits. The line emis- 
sion around 19 A is a mix of O VII He/3 emission and O VIII Lya emission. Lines 
from Fe-L transitions, in particular Fe XVII which has lines around 15 Aand 17 A are 
conspicuously absent. This may be surprising given the fact that SN 1006 is a likely 
Type la remnant, which produce ~ 0.7 Mq of Fe. However, for the low ionization age 
of SN 1006 Fe is expected to have a peak ionization state of Fe XIII, so no Fe XVII 
emission is to be expected. 
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4: The XMM-Newton-RGS spectrum of the northwestern knot. 
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FIGURE 5. The XMM-Newton PN and MOS spectra (left) and a detail of the Si emission as observed 
by the MOS detectors (right). The shifted line shows the position of the Si emission if most of the Si 
would have been in the He-like state. 



The CCD spectra were modeled using the SPEX^ [29] spectral code. The detailed 
results can be found in [19] and confirm the low ionization age of the plasma, net= 
2.4 X 10^ cm^-^s, whereas kTe= 1.5 ±0.2keV. There is a relative overabundance of Si, 
which suggests that the knot is an ejecta knot, and not a shocked pre-existing density 
enhancement in the interstellar medium. Interestingly, there is a weak enhancement of 
X-ray emission in front of the knot as well, that seems to emit slightly harder X-ray 
emission. It is not clear how these relate to each other. 

The low ionization age of the plasma can be deduced not only by globally fitting the 
spectrum, but also by measuring the centroid of each individual line, or by modeling 
the O VII triplet. For instance the EPIC PN and MOS spectra clearly indicate that the 
centroid of the Si emission is 1.80 ±0.01 keV instead of at 1.85 keV, which is the 
centroid expected if Si is in the He-like state (Fig. 5). Also the Ne and Mg line centroids 
are at lower energies than the He-like. Purely combining the information of the centroids, 
and assuming kTe= 1.5 keV one also finds an ionization age of Hg? ~ 2 x 10^ cm^^s 
(Fig. 6). The most important inner shell lines are present in the SPEX spectral code, but 
in general one should be aware that spectral codes are not always complete in this sense, 
as most atomic databases were initially compiled with CIE plasmas in mind. Moreover, 
the line parameters may also be less accurate, as less experimental data exists to measure 
the line properties. 



Non-equilibration of temperatures in the northwestern knot 

One of the main goals of the XMM-Newton observations was to measure the amount of 
non-equilibration of electron and ion temperatures, by measuring the thermal Doppler 
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FIGURE 6. Line energy versus charge for O, Ne, Mg, and Si. The crosses give the expected line 
energies for different ionization states (from Dr. M.F. Gu's FAC code) and the arrows approximately 
indicate the observed centroid energies. 



broadening of O VII and O VIII line emission. So far evidence for non-equilibration 
in supernova remnants comes mainly from measuring the Ha line width and broad to 
narrow line ratios. This has been done for several remnants including SN 1006 [30, 31]. 
The amount of non-equilibration has also been measured in the UV from C IV, N V 
and O VI lines by HUT [32] and very recently by FUSE [33]. The optical and UV 
results indicate that there is no equilibration at the shock front, although the Fi/SE' results 
suggests that the temperature of different species is not proportional to the species mass, 
as implied by eq. 1 and the HUT spectrum. 

The O VII Hea emission consists of the He triplet, plus a potential contribution from 
an O VI line very close to the O VII forbidden line emission. In order to measure the 
Doppler broadening models for the line ratios for a grid of n^t values were fitted, using 
parameters obtained from the flexible atomic code (FAC) of Dr. M. F. Gu [34]. The 
spectrum of the RGSl O VII spectrum and the best fit models are shown in Fig. 7. The 
best fit suggests a very significant broadening at a confidence level of 9.8a. Allowing 
full freedom of the line normalizations still indicates thermal broadening, at a still high 
confidence level of 6.5c7. The broadening itself is Oe = 3.4 ±0.4 eV. The published 
results are based on a response matrix that was corrected for the emissivity profile 
shown in Fig. 3. However, the fitting program XSPEC also contains the convolution 
model rgsxsrc made by Dr. A. Rasmussen, which uses an input image to determine the 
emissivity profile. Using rgsxsrc gives very similar results. 

The derived value for the broadening implies an oxygen temperature of kT = 530 ± 
150 keV, which looks incredibly high, but is in effect consistent with eq. 1 if the gas 
was shocked with a velocity of ~ 5000 kms^^ Of course the oxygen temperature is 
more than a factor 300 higher than the electron temperature, clearly indicating non- 
equilibration of temperatures. The measured O VII line broadening is more than reported 
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FIGURE 7. The OVII triplet from the northwestern not as observed by the RGSl. The dotted line 
shows the best fit spectrum with (sohd) and without including thermal Doppler broadening. 

for O VI UV emission from other regions of the northwestern rim [33]. There are 
several reasons why there could be a discrepancy. The knot may be an ejecta knot, 
as implied by the high Si abundance, and the heating may have been caused by the 
reverse shock, which has a different speed than the shock. Another possibility is that 
the plasma of the knot was shocked earlier than the regions observed in O VII, after all 
the X-ray measurements concern more highly ionized material, which is situated further 
downstream from the shock. 



What about the eastern rim? 

The eastern and western rims of SN 1006 are likely sites of efficient cosmic ray shock 
acceleration [26]. If the cosmic ray pressure dominates over the thermal pressure, eq. 1 is 
no longer valid. It will result in an average plasma temperature that is lower. Moreover, 
if the cosmic ray is very efficient and the cosmic ray power law index is equal to or 
flatter than -2, the energetic cosmic rays that escape the shock will drain the shock of 
its energy and the shock will become (cosmic) radiative, resulting in much higher shock 
compression ratios [35]. 

Given the interesting physics at the eastern and western rim, it is unfortunate that 
the RGS spectrum from the eastern rim does have a high spectral resolution, especially 
since no UV emission lines from the eastern rim have been detected by FUSE [33]. It 
is therefore still worthwhile to obtain as much information from the RGS spectrum as 
possible. I will therefore discuss here a preliminary analysis of RGS spectrum. 

The effective spectral resolution is mainly determined by the width of the shell, about 



Unfortunately in [33] the measured broadening in (7 is mistakenly taken to be a FWHM. 



5', resulting in an effective resolution for the first order spectrum of ~ 0.6 A. This is not 
sufficient to resolve the O VII triplet, which is also for the eastern rim the dominant line 
emission. However, there is still high spectral resolution information in the spectrum, 
which results from the sharp increase of the emission at the edge of the remnant. 

The O VII triplet emission is displayed in Fig. 8. It shows that a model that does 
only contain the triplet lines and an instrumental response that only takes into account 
the spatial emission profile does not give an adequate fit. The long wavelength side 
seems to fit reasonably well, but there is a dip in the data around 21.5 A, and the model 
underpredicts the emission around 21.3 A. Note that here the interior of remnant is 
dispersed in an opposite direction compared to the spectrum of the northwestern knot. So 
the edge of the remnant determines the shape of the long wavelength side. This suggests 
that discrepancy is caused not so much by line broadening, but by line shifts. The plasma 
motions at the edge of the remnant, will be mostly in the plane of the sky, whereas going 
to the interior of the remnant radial motions will give rise to Doppler shifts. As we are 
observing both the front and backside of the shell, both blueshifts and redshifts should be 
taken into account. The redshifts are not directly noticeable (they blend in with the rest 
of the triplet emission), but the blueshifts can be seen as the mismatch between model 
and spectrum around 22.3 A. 

I have attempted to fit a simple model for the O VII triplet, taking into account the 
increase in blueshifts and redshifts moving interior to the remnant (corresponding to 
shorter wavelength). As it is a mixing of spectral and spatial information, the method 
consisted in making a grid of response matrices in which the spatial and redshift and 
blueshift were incorporated. The input to the model consisted of the projected plasma 
velocity, which increased toward the interior, and a parameter that determined how 
much of the emission was coming from the near and far side of the remnant. This is 
a simplification, as velocity gradients are presents and the morphology of the shell may 
be complicated. It was nevertheless possible to obtain acceptable fits to the O VII triplet 
emission (Fig. 8). 

It turns out there are two good fitting models. One with no, or only a small additional 
(thermal) broadening (Ge = 0.4 ±0.3 eV), and a shell with velocity 6500 ± 830 km s~^ 
(68% confidence range. Fig. 8). The other solution is for a much higher shell velocity 
of > 9000 kms~\ and quite a substantial additional broadening of (Te = 3 ± 1 eV. In 
both cases the best solution suggests that about 70% emission comes from material 
moving toward us (> 55% at a 3a confidence). The solution with a large broadening 
gives a better fit (Fig. 8), but is less plausible. The reason is that a bulk velocity of 
> 9000 km s~^ would imply that the shock heated plasma is still in free expansion, 
whereas a solution with v ~ 6500 km s~^ fits in much better with radio proper motions 
[36]. The average shock velocity of SN 1006, given its age and distance is < v >= 
9375 km s^^ The radio expansion measurements suggest that = (0.48 ±0. 13) < v >= 
(4500 ± 1220) km s~^ Both solutions exclude the possibility that the plasma velocity is 
similar to what has been measured for the northeastern Ha filaments Vv = 2800 km s~^ 
[30]. This is consistent with the smaller radius in the northwest and the finding of an 
increased density as indicated by 21cm line measurements [37]. 

The shape of the O Vn triplet suggests that more material is blueshifted than red- 
shifted. This asymmetry is reminiscent of the difference in forward shock and reverse 
shock velocity found in absorption toward a UV star projected on the center of SN 1006 
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FIGURE 8. The O VII triplet from the eastern rim as observed by RGS 1 . Top left: The best fit model 
with only taking into account the spatial profile. Top right: the likelihood statistic as a function of shell 
velocity, for a front back asymmetry of 70%. The arrows indicate other velocity measurements from 
respectively the northwest, radio and the average velocity. Bottom: The left panel shows the best fit 
solution for a shell velocity of 6500 km s^' with no broadening. The right panel shows a solution with 
a velocity of 9500 km s^^ including substantial thermal broadening. 



[38]. The reverse shock appears to have penetrated more into the innerejecta layers of 
the front side than on the backside. Such an asymmetry is likely to be the result of a 
higher density at the front side, which will also result in more shocked heated plasma at 
the front side. 



CONCLUDING REMARKS 

SN 1006 is from many points of view an interesting supernova remnants. It is one of 
the few historical remnants, and it was the first remnant for which it was established 
that the particle acceleration by the shock was efficient enough to produce X-ray syn- 
chrotron emission. The fact that SN 1006 expands into a low density medium makes 
it also an ideal laboratory to study the effects of non-equilibrium ionization and non- 
equilibration of temperatures, both of which depend on the value of n^t, which is around 
2 X 10^ cm^s for SN 1006, one of the lowest values found in a supernova remnant. 



The measurements presented here show the diagnostic power of high resolution spec- 
troscopy when it comes to characterizing the ionization state, temperature and velocities 
of hot plasma. The results on a remnant like SN 1006 raises the expectation for the re- 
sults from micro-calorimeters carried hyAstro-E2, and Constellation-X/Xeus, which will 
make it much more straightforward to measure Doppler motions from extended sources. 

I thank the organizers of the XDAP 2004 workshop for their invitation and financial 
support. 
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